Entropy, frustration and large thermopower of doped Mott insulators on the fee lattice 
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Electronic frustration and strong correlations may lead to large Seebeck coefficients. To understand 
this physics on general grounds, we compute the thermopower of the one-band Hubbard model on 
the 3-dimensional fee lattice over the whole range of fillings for intermediate and large interaction 
strength. Dynamical mean-field theory shows that when the density approaches half-filling, the fee 
lattice at strong coupling exhibits a large low temperature Seebeck coefficient S. The largest effect 
occurs as one approaches n — 1 from dopings where electronic frustration is maximized. The high- 
frequency limit of the thermopower and the Kelvin limit are both used to provide physical insight 
as well as practical tools to estimate the thermopower. The high-frequency limit gives a reliable 
estimate of the DC limit at low temperature when the metal becomes coherent. By contrast, the 
Kelvin approach is useful in the strongly interacting case at high temperature when transport is 
incoherent. The latter result shows that in doped Mott insulators at high temperature and strong 
coupling the thermopower can be understood on entropic grounds. 
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I. INTRODUCTION 

The search for better thermoelectric materials has 
focussed on strongly correlated materials in recent 
years, going beyond traditional semiconductor systems. 
Narrow-band correlated oxides in quasi two dimensio- 
nal systems such as sodium cobaltate Na x CoC>2 1 show 
large thermoelectric power (or Seebeck coefficient S) S ~ 
100fiV/K, at ~ 100K. Mott Hubbard physics is invoked 
to understand this extraordinary scale of S in Na x Co02 
and Sr x La\^ x Ti0^ 2 , often in the limit of zero bandwidth 
or atomic limit t — > 0, 2 ~ 4 . More recently, FeSb2 (a Kondo 
insulator or a correlated band insulator) has drawn much 
attention. 5-7 This is a three-dimensional compound with 
a colossal Seebeck coefficient S ~ 45,000 /J.V/K that is 
of great interest for thermoelectric cooling at cryogenic 
temperatures. 

For more traditional classes of thermoelectric mate- 
rials, such as the semiconductor Bi^Te-i, the Bloch- 
Boltzmann transport theory points to a way to optimize 
the band structure for thermoelectric applications 8 . In 
the case of strongly correlated materials, the situation is 
much more complicated due to the fragility of the quasi- 
particles. These appear with strongly renormalized physi- 
cal properties only in the so-called coherent regime that 
may develop at very low temperatures in the presence 
of correlations. Strongly correlated systems favour large 
thermopower then because the large effective masses cor- 
respond to small Fermi temperatures. The temperature 
scale for maximum thermopower is set by fractions of 
hopping t. Estimates for real materials give for V2O3 9 
t « 58CLK", for FeSi, 10 , t w 152 Jf, and t « lOOif for 
Na^CoC^ 11 . Hence the peaks that are found for T < t 
as well as the strong temperature dependence up to about 
T ~ t are in a physically relevant temperature range not 
far from room temperature. Room temperature may even 
fall in the incoherent regime where Fermi liquid theory 
does not apply and S can be large. 



Early studies of the effect of correlations on the ther- 
mopower focused on the Hubbard model on unfrustrated 
lattices 12-16 . On these unfrustrated lattices, it is unli- 
kely that large values of thermopower S can be achie- 
ved 17 because when the temperature T is of the order 
of the coherence temperature T coh , the thermopower has 
contributions from the quasiparticle peak but also from 
the adjacent Hubbard band and these contributions have 
opposite sign. Study of the effect of frustration on the 
cubic lattice with next-nearest-neighbor hopping at half- 
filling for interaction strengths below the Mott transi- 
tion 18 show that the peaks that are obtained in S are 
similar to those in the electronic specific heat and are 
associated with the thermal destruction of the quasipar- 
ticles. More recently, realistic calculations 17,19 of ther- 
mopower for experimentally interesting compounds have 
shown increase in thermopower due to correlations 7 ' 20 ' 21 . 
The role of frustration has been emphasized in particu- 
lar for the cobaltate Na x Co02 using the t — J model on 
a 2d finite triangular lattice 11 . Realistic calculations for 
that compound 22 ' 23 demonstrate that frustration and in- 
teractions combined also with disorder can lead to large 
enhancements. 

Here we focus on a generic model system with no di- 
sorder. More specifically, we study the fee lattice be- 
cause of its prevalence in materials and because its 
electronic structure is frustrated already for nearest- 
neighbor hopping. Shastry and coworkers 24-27 have ar- 
gued that in correlated materials, electronic frustration 
can amplify the Seebeck coefficient considerably beyond 
Bloch-Boltzmann (band) theory expectations. The fee 
lattice has optimal thermopower, according to the Bloch- 
Boltzmann theory 8 , when the chemical potential coin- 
cides with the sharp peak in the density of states. Here 
we show, instead, that in the presence of interactions, 
large thermoelectric effects appear at low temperatures 
close to half-filling and for dopings that correspond to 
large electronic frustration, provided the system is suffi- 
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ciently correlated that the half-filled band is insulating. 
We show that insights on the temperature and doping de- 
pendence of our results may be obtained from methods 
developed recently by Shastry and coworkers 24-27 . In 
these approaches, one considers either the high-frequency 
limit or the ultra- low frequency (Kelvin) limit. In the 
latter case, the thermopower is related to the density 
dependence of the entropy at fixed temperature. We 
find that such entropic considerations beyond the Bloch- 
Boltzmann theory can control the thermopower in the 
incoherent regime and that in the coherent regime the 
high-frequency limit is often reliable. In other words, 
while approximate, these approaches also give correct 
order-of-magnitude estimates and bypass the difficulties 
facing Kubo formula calculations. These methods can be 
straightforwardly implemented in existing codes for rea- 
listic calculations. 28 

In the next section we describe the model and the Dy- 
namical mean-field theory (DMFT) approach along with 
the three formulas used for thermopower : Kubo, high- 
frequency and Kelvin. The results are then discussed for 
two values of interaction strengths, half bandwidth and 
two times the bandwidth, in the context of the Kubo for- 
mula. The next to last section makes comparisons with 
the high-frequency and Kelvin approaches. We end with a 
summary. Appendix A shows that in the strong-coupling 
limit the thermopower changes sign at half-filling n = 1 
instead of n = 0.8 in the Bloch-Boltzmann theory. Ap- 
pendix B presents benchmarks, in the worse case, of the 
impurity solver we used for the Kubo calculations. 



II. MODEL AND METHODS 
A. DMFT 

Dynamical mean-field theory (DMFT) has proven its 
accuracy for the Hubbard model and its predictive po- 
wer for real materials 29 ' 30 that have a clearly three- 
dimensional structure, despite the fact that the me- 
thod is exact only in infinite dimension. Here we ap- 
ply it to the one-band Hubbard model on the fee lat- 
tice where the single-particle dispersion is given by £k = 
— 4i[cos(fc :c ) cos(ky) + cos(k x ) cos(/c z ) + cos(k y ) cos(fc z )]. 
We take nearest-neighbor hopping t positive. 

We consider only the paramagnetic phase. To solve the 
impurity problem of DMFT, we use two different me- 
thods. 1) The numerically exact continuous time quan- 
tum Monte Carlo method (CTQMC) 31 , a finite tempera- 
ture approach that relies on the Monte Carlo summation 
of all diagrams obtained from the expansion of the par- 
tition function in powers of the hybridization A. This 
method does not have errors associated with time discre- 
tization or bath paramctrization and is therefore exact 
within statistical errors, but computationally expensive. 
When analytical continuation is an issue for this method, 
we rely on : 2) Iterated Perturbation Theory (IPT), an 
approximation method based on an interpolation from 



2 nd order perturbation theory for the Anderson impurity 
problem 32 . The interpolation preserves the correct high- 
frequency limit for the self-energy and is exact in both 
the non-interacting and the atomic limits. For U jt in the 
intermediate-coupling regime, the condition of Rcfs. [33] 
and [18] to fix the bath occupation number leads to ac- 
curate results. However, at finite temperature and large 
U, we need to use the version of IPT that was developed 
and carefully benchmarked in Rcf. [34] . In this approach, 
IPT-D, one needs double occupancy. At very large U, 
the simplest estimate for double occupancy suffices. We 
checked for many observables that CTQMC and IPT in 
this particular implementation give similar temperature 
dependent properties. 



B. Thermopower 

A general Kubo formula for the thermopower may be 
written down using linear response theory 24 as 



S(q x ,w) = 



(1) 



where the susceptibility of any two operators A, B is gi- 
ven by X a,b(!») = *Jp°° dt f"-° +t <M*), B(0)]> and 
where p, K = H — pN and J x arc the charge density, 
the (grand canonical) Hamiltonian and the current ope- 
rator respectively at finite wave vectors. The fast limit 
where q — > first and then u — > 0, is the relevant li- 
mit, that we will call Sxuboi for DC transport. 35 In the 
DMFT limit the vertex corrections vanish for one band 
models 36 ' 37 but one still needs analytical continuation of 
single-particle Green's functions since the methods of so- 
lution for DMFT are formulated in Matsubara frequency. 
That may be done to obtain meaningful results in the 
simplest cases 18 . However, the most accurate modern me- 
thod to solve the DMFT equations, namely CTQMC, 
requires analytical continuation of numerical data that 
contains statistical uncertainties. The fact that the ker- 
nel for Xj/ q \ K{- q!C )(. u ) 1S not positive definite renders 
this in practice extremely difficult, even for maximum en- 
tropy methods that are usually called for with CTQMC 
results. That is why to compute SKubo we performed the 
analytical continuation on the self-energy obtained by 
IPT-D using the Pade N points algorithm that can be 
found in Ref. [38] . Sxubo is then calculated using the real 
frequency equations 18 . In Ref [34], we showed that the 
results obtained for the resistivity were at least qualita- 
tively valid. 

Two new approximate methods to compute the 
thermopower 24,27 without analytical continuation, na- 
mely the superfast limit S* and the Kelvin formula 
S Kelvin, are explained below. They have the advantage 
that they do not require analytical continuation. In addi- 
tion, they have physical content that will help us identify 
the origin of enhanced thermopower. 
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Consider first the superfast limit of Eq.(l), S* = 
lim w>W( , j^—to S(q, uj), where u> c is the largest characteris- 
tic frequency in the problem. In the case of the Hubbard 
model, we would identify lu c = max{W, [/}, where W is 
the bandwidth. In terms of the thermoelectric operator 
§ xx and the stress tensor t xx that appears in the f-sum 
rule , the quantity S* is equal ^pss\ ■ Since the Seebeck 
coefficient is a ratio of transport coefficients where the 
scattering rate cancels in the relaxation time approxima- 
tion, it is conceivable that it may have a weak dependence 
on frequency. It turns out that S* captures much of the 
many-body content of the exact SKubo when the charac- 
teristic u) c is not too large. 24 For the Hubbard model, 
the commutators can be evaluated exactly and the cor- 
relation functions evaluated solely in terms of the single- 
particle Green's function 28,39 
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Figure 1: (Color online) S in /jV/K for U = as a function 
of T/t for different values of density : n = 0.2 (yellow (<)) 
0.4 (black (o)), 0.6 (blue (*)), 0.8 (red (□)), 1.2 (green (0)), 
1.4 (cyan (V)), 1.6 (magenta (A)) and 1.8 (brown (*)). 
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with q e the (negative) charge of the electron. While one 
can show that lim-r-^o^ xx ) — when U — 0, for interac- 
ting systems it was found in numerical calculations 25 ' 26 ' 28 
and from Fermi liquid theory that ]imr_+o{<f> xx ) an d 
hence, S* diverges as T vanishes. Indeed, at large fre- 
quency there are reactive parts to the energy transport 
that should not be considered part of the thermal current. 
We thus eliminate the constant term limT->o(3? xa: ) = 
(<j> xx )o by fitting the low temperature part of (§ xx ) to 
(<j> xx )o + bT 2 , a functional form derived from the Som- 
merfeld expansion that empirically remains valid in the 
presence of strong interactions. We then define the ther- 
mopower in the super fast limit by 



S* 



1 



T 



(t xx ) 



(2) 



Not subtracting {§ xx )o leads to large deviations between 
S** (now only S*) and SKubo at low T for all dopings, 
contrary to what we find. 28 

The second approximate method considers the slow li- 
mit of the general formula for thermopower, Eq. (1). As 
noted earlier 24 ' 27 , this leads to a finite and interesting 
answer that is termed the Kelvin formula 



SKeivin = lim S(q,u>) = — 

9->0 q e V on 



q e 



da 
df 



T,V y e \"-"- / V,n 

(3) 

where s is the entropy density, n the number density, V 
the volume, fj, the chemical potential, ks Boltzmann's 
constant and where the last form follows from a Maxwell 
relation. The formulas for S** and SKeivin approach 
SKubo from two different limits. 
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Figure 2: (Color online) SKubo, SKeivin and S** as a func- 
tion of T/t as calculated with IPT for Kubo and CTQMC 
for the Kelvin and high-frequency estimates for two values of 
interaction U = 8t (a), (b), (c) and U = 32t (d), (e), (f) for 
different values of density : n = 0.4 (black (o)), 0.6 (blue (*), 
0.8 (red (□)), 0.84 (black - - (<)), 0.88 (blue - - *), 0.92 (red 
- - (x))), 1.08 (green - - (>)) 1.2 (green (0)), 1.4 (cyan (V)) 
and 1.6 (magenta (a)). 
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III. RESULTS 
A. Band limit 

The Seebeck coefficient in the band limit (U = 0) and 
constant relaxation-time approximation is displayed in 
Fig. 1 along with the non-interacting density of states 
as an inset. The bare bandwidth is 16i. The horizontal 
triangles D> on this plot indicate the temperature Tg, be- 
low which the leading term of the Sommerfeld expansion 
for the U = case is 90% of the full answer. The sign 
and magnitude of the results in Fig. 1 depend mostly 
on the particle-hole asymmetry (X(oj) ^ X(—u>)) of the 
energy-dependent transport function 

*H = E(|;)V- £k ) (4) 

computed in Ref. [34]. At the qualitative level, the be- 
havior of the thermopower is easy to understand. Below 
Tg it is linear, the behavior expected from Fermi liquid 
theory. The low temperature thermopower changes sign 
from negative to positive when the Fermi surface turns 
from electron-like to hole-like at the van Hove singula- 
rity in the density of states, located around n = 0.8. In 
the large temperature limit, the thermopower saturates 
to a value of the order of the fundamental unit of ther- 
mopower, fcs/e = 86[iV/K. As can be seen from the 
sharp peak in the density of states, the absolute value 
of S is maximum in the large T limit and for an almost 
filled band where the largest deviations from particle-hole 
symmetry occur. This asymmetry comes from the fact 
that the electronic spectrum for large filling is frustra- 
ted : the extrema of the band dispersion for excitations 
arc not at high-symmetry points in the Brillouin zone. 
This may be understood simply from the electronic spec- 
trum of a single triangle. Any site forms, with any pair of 
its nearest-neighbors on the FCC lattice, an equilateral 
triangle. Note that frustration leads to low values of Tg. 

B. Weak and strong coupling limits 

Armed with this understanding, we turn to the two in- 
teracting cases displayed in Fig. 2. The value U = 8t, 
namely half the bandwidth in Figs. 2(a),(b),(c), is in 
the intermediate coupling range whereas U = 32i, twice 
the bandwidth, in Figs. 2 (d),(c),(f) is representative of 
strong coupling. These values of U are respectively below 
and above the Mott transition at half-filling. For each va- 
lue of U we evaluate S in three different ways : SKubo, 
Sneivin and S** . 

Let us move slowly through the many different results 
displayed in Figs. 2. In these figures we do not display 
the results far away from half-filling, namely n = 0.2 and 
n = 1.8, since they are essentially the same as in the band 
limit Fig. 1. Indeed, at low density the carriers (electrons 
at n = 0.2 and holes at n = 1.8) can avoid each other 



and the interaction is not important. Raising the carrier 
densities, we see that for n — 0.4 (black) and n = 1.6 
(magenta) the results still do not depend much on U . In 
addition, the results for SKubo and S** are very close to 
each other. 

As we move closer to half-filling, the thermopower be- 
comes more sensitive to interactions. Let us then first 
focus on SKubo Fig. 2(a) at weak to intermediate cou- 
pling, U — 8t. At densities n — 0.6 (blue) n = 0.8 (red) 
and n = 1.2, (green) n = 1.4 (cyan) the main qualitative 
changes are a sign change as a function of T for n = 0.6 
(blue) n = 0.8 (red) and the appearance of a small low 
temperature peak for the four fillings. However, the ab- 
solute values do not exceed the largest absolute values 
that can be reached in the band limit by doping. 

At strong coupling, U — 32i, there are larger devia- 
tions from the U = results. The system is an insulator 
at n = 1. Upon doping, it develops a large Seebeck coef- 
ficient with non-trivial T dependence even below Tg defi- 
ned in Fig.l. In Fig. 2(d), at the fillings discussed above, 
n = 0.6 (blue) n = 0.8 (red) and n = 1.2, (green) n = 1.4 
(cyan), the peaks in SKubo are more pronounced at large 
U and there is a sign change at finite temperature even 
for n = 1.2. The maximum magnitude of thermopower 
is obtained at low temperature : for example for n = 0.8 
(red), S « -15/j.V/K at T = 0.3t, or for n = 1.2 (green), 
S w 57/j.V/K at T = 0.24i, or for n = 1.08 (dashed 
green), S w 71.23^U/A' at T = O.Olt. This represents at 
least an increase by a factor 3-4 to more than an order of 
magnitude with respect to the non-interacting case. The 
appearance of pronounced peaks at low temperature that 
are absent in the non-interacting limit reflects, as pointed 
out in Ref. [18], the fact that Fermi liquid coherence is 
recovered only at very low temperature in the presence 
of strong correlations. 29 The peaks for S are similar to 
those of the specific heat (see Fig. 9 of Ref. [34]) and de- 
fine T coh , as noted before 18 . Since T coh in IPT-D and 
CTQMC is similar 34 the positions of the peaks in SKubo 
should be correct. Frustration is a factor that contributes 
to make the Fermi liquid coherence T low. It also plays a 
role in the t — J model on the triangular lattice. 11 . The 
results also confirm the predictions of Shastry using S** 
concerning the dependence of the thermopower on the 
sign of t on a frustrated lattice. 

The results very near half-filling, shown by dashed lines 
in Fig. 2(d) are especially interesting : a) Large ther- 
mopower develops at low temperature above half-filling 
where there is frustration, b) The low temperature slope 
changes sign near n = 1 by contrast with the smaller U 
results where the sign change occurs near n — 0.8. We 
show in Appendix A that the change of sign of SKubo 
at low temperature near n = 1 can be confirmed di- 
rectly from the CTQMC results. From a Fermi liquid 
point of view this phenomenon reflects the importance of 
particle-hole asymmetric terms ujT 2 and u 3 in the scat- 
tering rate. 17 
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IV. OTHER MEASURES OF THERMOPOWER : 

Here we compare S**(T) and SKeivin{T) respecti- 
vely in Figs. 2(b), (e) and (c),(f) with the corresponding 
SKubo(T) . The results for S**(T) and SxeiviniT), obtai- 
ned directly with CTQMC without analytical continua- 
tion, provide both physical understanding and an alter- 
nate more straightforward way to estimate the thermo- 
power. For S**(T), results for both U = 8t and U = 32t 
obtained using IPT give basically the same results as 
those obtained from CTQMC. Hence they can be consi- 
dered reliable. For SKeivin(T) at large U = 32t, we ob- 
tained the same qualitative behavior of the thermopower 
but with smaller coherence temperatures for CTQMC, as 
shown in Appendix B. Above the coherence temperature, 
even the quantitative differences disappear. 

We focus first on the low T linear regime, which cor- 
responds most of the time to operating conditions. In 
the spirit of the Boltzmann calculation for U = in the 
constant relaxation-time approximation, we have a phe- 
nomenological form for S at low T involving two different 



Fermi temperature T. 



and r} 2) 



S ~ Sgn(S) 



k B it 2 T 
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(i) 



(5) 



and 



Depending on the band structure and filling, Tj 

(2) 

Tj can be very different from each other even at U = 0. 
We wrote explicitly the sign of the thermopower to be 
able to define a positive T^K Fig. 3 shows Sgn(S)/T^ 
as calculated from the Boltzmann non interacting S (so- 
lid black line), and for Sxubo (blue squares). The latter 
is our reference value for finite U. For comparisons, we 
plot S** (red circles), and S Kelvin (brown diamonds) as 
well as Sfl (green star), the value obtained by assuming 
that we have a Fermi liquid with a pure uj 2 dependence 
without deviations from particle-hole symmetry, i.e. the 
first term of S Fermi in Eq.(A4). 17 The left panel is for 
U = 8t (a) and the right panel for U = 32t (b). This 
illustrates two of the important results of this paper i.e. 
(1) for weak to intermediate coupling in Fig. 3 (a), S** 
is a good indicator of the behavior of S. (2) The failing 
of Spceivin in that region suggests that the common view 
of the thermopower as the entropy per particle is not 
to be taken seriously when the system is coherent. Still, 
Sxeivin for large U, Fig. 3(b), predicts the correct sign for 
the low T thermopower. Fig. 3(b) also shows that S**(T) 
gives information about the increase of S close to n = 1 
even if its wrong just below half-filling where the quasi- 
particles weights become very small while the scattering 
rates become large 34 and particle-hole asymmetric. By 
comparing Sfl (green star) first term with Sxubo (blue 
squares), we finally see that at strong coupling the effects 
of the non-quadratic terms of the self-energy on the low 
T thermopower are important mostly for n < 1. 

The fact that S** is less reliable for the Hubbard model 
at strong coupling, where interaction matters, contrary 
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Figure 3: (Color online) Sgn{S)/Tf ] for U = 8t (a) and 
U — 32t (b) as a function of the density. The Boltzmann non- 
interacting value is given by the black solid line, SKubo (blue 
(□)), SFermi with only the first term (green (*)), Sneivin 
(brown (0)) and S" (red (o)). 



to the t — J model, 11 , was expected from the fact that 
large U makes cj c , the largest characteristic frequency in 
the problem, too far from the u> = limit. 

Overall then, one finds from the results in Fig. 2 that 
for weak to intermediate coupling (U = 8i), the frequency 
dependence of the scattering rate is not so important and 
S** gives a good estimate of the thermopower. On the 
other hand, for strong coupling (U = 32t) the situation 
is more complicated. At low T, where single-site DMFT 
gives a Fermi liquid regime, S** seems to give a good 
estimate when the carrier density is not too close to half- 
filling where interactions effects become important. By 
contrast, when Fermi liquid behavior is present S Kelvin 
is not reliable at low temperature and strong coupling. 
At high T however, SKeivin gives qualitatively correct 
results, 40 accounting for the fact that the temperatures 
where the thermopower changes sign are off. This de- 
monstrates, through Eq. (3), that in the strong coupling 
incoherent regime, entropy drives the thermopower. The 
relationship between entropy and thermopower has been 
discussed in the linear low temperature regime and mul- 
tiband systems experimentally 41-43 and theoretically 44 . 

Let us come back to the high temperature regime, 
U > T > t. In this limit, the so-called Mott-Heikes 
formula also suggests a relationship between the en- 
tropy per particle and thermopower. Heike's formula, 
S —t jJb/T = (ds/dn)E,v is derived by applying the ato- 
mic limit to the case of large U. The thermopower in this 
case depends only on density since \xjT is T independent. 
It is given by 45 



Smh = 
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, 2-n ) 



if < n < 1, 
if I < n < 2. 



(6) 



One striking prediction of these formula is a change of 
sign of the thermopower at 1/3 doping from half-filling. 
None of our results in Fig. 2 are in the regime where \xjT 
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is T independent 34 . It is quite remarkable then that for 
temperatures between T = t and T = 2t we nevertheless 
have more than one sign change as a function of density 
for a fixed T. 

Even though S Kelvin = Smh in the regime where /j/T 
is T independent, when T is not high enough for this 
condition to be realized, as in our case, SKeivin makes 
a stronger case than Heikes formula. Indeed, SKeivin = 
(ds/dn) T y relies on the correctness of interchanging the 
order of two limits, which can become justified in the in- 
coherent regime 46 . It was checked explicitly in Ref. [40] 
that Kelvin's formula is better than a modified Heike's 
formula, (/x(T) — /x(0))/T. To illustrate what happens in 
our case, we show in Fig. 4 the thermopower as a function 
of density for T = 2t and U = 32t for SKubo, S Kelvin, S** 
and compare with Smh in Eq. (6). For densities above 
about n = 0.8, our Kubo results match quite well with 
Heike's. For large doping, S** works better than SKeivm 
or Smh, which is another indication that large interac- 
tion effects are necessary for S Kelvin or S mh to be good 
approximations. Despite their simplicity, we see that the 
latter two formula capture the essence of the thermopo- 
wer in the high temperature regime. Note that the quan- 
titative results are a bit off, especially for n < 0.8 where 
the absence of electronic frustration seems to leave more 
room for coherence. At this particular temperature, Hei- 
ke's formula seems to be somewhat better overall than 
Kelvin's. However, in the presence of electronic frustra- 
tion at n > 1, Heike's formula shows a small systematic 
deviation from Kubo's that does not appear for S Kelvin- 

Chakraborty et al. ' studied the case of a square lat- 
tice with nearest-neighbor hopping at large U = 20t for 
T = 0.5t in the incoherent regime. U = 40t and T = t 
would be more comparable to our case given that the 
bandwidth of the fee lattice is twice as large as that of 
the square lattice for the same t. They solved the pro- 
blem using CDMFT with NCA as an impurity solver. 
They did the calculation for densities between n = 0.6 
and n = 1 and found, like we do, that Kubo's result is 
quite close to Heike's. A more thorough investigation of 
the high-temperature regime is called for. 
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Figure 4: (Color online) Si< u bo (black -), S Kelvin (blue - - 
), Smh (red -.) and S** (magenta o CTQMC) in fiV/K for 
U — 32t as a function of density for T = 2t. In order to 
calculate Sneivin we used IPT results. As shown in Fig. 5, 
in this temperature range, IPT-D gives essentially the same 
fi(T) as CTQMC. 



tronic frustration. We leave open the question of whether 
there is an optimal doping or just a sudden jump from 
very large thermopower at infinitesimal doping to zero at 
half-filling. For doped Mott insulators in the incoherent 
large-temperature limit, the Kelvin formula gives a rea- 
sonable estimate of the thermopower without the need 
for analytical continuation. This clearly shows that large 
thermopower in this case can be explained on entropic 
grounds over the whole doping range at large T. 

It would be of great interest to verify how SKeivin be- 
haves at low T and large U when non- Fermi liquid beha- 
vior is observable. For this, one needs to consider lower- 
dimensional systems that require going beyond single-site 
DMFT. Also, addition of disorder trough CPA 48 or sta- 
tistical DMFT 49 would be of interest to verify how S** 
and SKeivin inform us about the interplay between cor- 
relations and disorder. 



V. SUMMARY 

We have shown that from weak to intermediate cou- 
pling, the thermopower of the FCC lattice is not in- 
fluenced very much by interactions. The quantity S**, 
that neglects the frequency dependence of the scattering 
rate, can be a reliable estimator for the thermopower. By 
contrast, at strong coupling, namely when at n = 1 the 
system is a Mott insulator, sharp peaks in the T depen- 
dence of Seebeck coefficient appear at low temperature 
and, for the range of dopings considered, absolute values 
comparable to ks / \ q e | , close to hundreds of (J.V/K, can be 
obtained. We found a sharp increase of the thermopower 
at a given T when one approaches n = 1. The increase is 
particularly sharp for n > 1 where there is strong elec- 
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Annexe A: Sign change of the thermopower around 

n = 1 for U = 32t 



Since the calculation of SKubo rests on IPT-Z), we need 
to verify that the change of sign of the low temperature 
slope around n = 1 is not an artefact of the method. Wc 
present three arguments. First, we show that it is possible 
to obtain the change of sign purely from CTQMC results. 
Although analytical continuation for CTQMC is nearly 
impossible for such a large U, there is a way to check 
the sign of the thermopower directly from the behavior 
in Matsubara frequency. Following Ref. [50] we define the 
function 



lij (fin) 



r — 

J-oo 71 



duj n n 



L 13 (uj). 



(Al) 

In Eq. (Al), i and j can take the value 1 or 2. \u 1S the 
charge-charge, X12 is the charge-heat and X22 is the heat- 
heat susceptibility. The sign of the thermopower is thus 
given by the sign of Li 2 (0). In terms of the 7^ function, 
we can show that 



M°) = a lim o 7 <j (0„). 



(A2) 



If we replace the discrete variable fl n by a continuous 
variable z, Eq. (A2) would correspond to 



lim iij{9. n ) -> lim 7y(2) 

f2„— S-0 2—5-0 



-Xij(z) 



(A3) 



Thus the DC value would be given by the slope at the 
origin. The problem is to find the correct continuous func- 
tion since the value of the derivative will depend strongly 
on the function. But for the sign there is no such problem. 
Direct inspection of the fi„ dependence indicates what is 
the sign of the derivative and thus the sign of the ther- 
mopower. To verify again this result, we also calculated 
the susceptibility X12 by the same method we used for 
X11 in Ref. [34] but using results obtained via CTQMC. 
We find that, indeed, at U = 32i, the low temperature 
thermopower changes sign around half-filling rather than 
around the small U result n — 0.80. Finally, we can also 
use just a few Matsubara frequencies to obtain the self- 
energy near zero frequency from Pade approximants and 
use it to calculate £12(0). The sign should be correct even 
though the actual value cannot be trusted. Once again 
we obtain that the sign change occurs around half-filling 
if we use 5, 10 or 20 Matsubara frequencies. We are thus 
confident that the IPT-D results for SKubo are at least 
qualitatively significant. 



The sign change and its relation to the underlying 
Fermi liquid predicted by DMFT can also be understood 
analytically. Starting from an u 2 self-energy, Refs. [12, 
17], demonstrated that the thermopower only depends on 
interaction through Z, the quasiparticle weight. Hence, 
at first sight the sign seems determined solely by the non- 
interacting problem, namely by the sign of the derivative 
of the transport function Eq. 4. Howewer, when one consi- 
ders the particle-hole asymmetry of the quasiparticle life- 
time by including wT 2 and w 3 terms in the self-energy 17 , 
the thermopower is given by 



k B k B T 



'Fermi 



0(A) El 



a 2 E 2 



(A4) 



where /2 = fj, — Re{X(0)}, the E's are universal constants 
while ai, a 2 and 70 are fitting parameters for the different 
powers of the self-energy and thus have an interaction de- 
pendence. Recent strong coupling theories 51 provide an 
analytical understanding of the odd frequency contribu- 
tions to the quasiparticle lifetimes, which become pro- 
minent as the particle density increases towards half 
filling. These provide dynamical asymmetry and signifi- 
cantly influence the ARPES line shapes. From Eq. (A4), 
we see that the asymmetry in the lifetime creates a term 
that might compete for the sign and magnitude of the 
thermopower. This is what is happening in our case. Be- 
low half filling but above n = 0.8, the band results would 
be positive but the interacting part leads to a strong 
particle-hole asymmetry that competes and keeps SKubo 
negative. This compensation leads to a small thermopo- 
wer between n = 0.8 and n = 1. However, above half- 
filling, either both contributions are of the same sign 
and enable an increase in the thermopower or the second 
contribution is negligible. The asymmetry in the scatte- 
ring rate also plays an important role in the temperature 
dependence of the thermopower, as shown in Ref. [40]. 



Annexe B: S Kelvin at U = 32t 

The only quantity that does not require analytical 
continuation and depends on the method of calculation, 
i.e. IPT-D vs CTQMC, is SKeivin at strong coupling. The 
comparisons are made in Fig. 5. Lines with the same co- 
lor are obtained for he same filling. There are dots when 
they are obtained from CTQMC. The main quantitative 
difference is that the coherence temperatures (peaks) are 
smaller with CTQMC and closer to the ones predicted 
by SKubo and C„ even if those two quantities were ob- 
tained with IPT-D. This quantitative difference does not 
prevent qualitative comparisons. For T larger than the 
peaks, then even the quantitative differences disappear. 
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